OPEN 3 ACCESS Freely available online 



•0-PLOS I ONE 



Neurospora crassa Protein Arginine Methyl Transferases 
Are Involved in Growth and Development and Interact 
with the NDR Kinase COT1 

Daria Feldman, Carmit Ziv, Rena Gorovits, Michal Efrat, Oded Yarden* 

Department of Plant Pathology and Microbiology, The R.H. Smith Faculty of Agriculture, Food and Environment, The Hebrew University of Jerusalem, Rehovot, Israel 



Abstract 

The protein arginine methyltransferaseas (PRMTs) family is conserved from yeast to human, and regulates stability, 
localization and activity of proteins. We have characterized deletion strains corresponding to genes encoding for PRMT1/3/5 
(designated amt-1, amt-3 and skb-1, respectively) in Neurospora crassa. Deletion of PRMT-encoding genes conferred altered 
Arg-methylated protein profiles, as determined immunologically. Aamt-1 exhibited reduced hyphal elongation rates (70% of 
wild type) and increased susceptibility to the ergosterol biosynthesis inhibitor voriconazole. In Aamt-3, distances between 
branches were significantly longer than the wild type, suggesting this gene is required for proper regulation of hyphal 
branching. Deletion of skb-1 resulted in hyper conidiation (2-fold of the wild type) and increased tolerance to the chitin 
synthase inhibitor polyoxin D. Inactivation of two Type I PRMTs {amt-1 and amt-3) conferred changes in both asymmetric as 
well as symmetric protein methylation profiles, suggesting either common substrates and/or cross-regulation of different 
PRMTs. The PRMTs in N. crassa apparently share cellular pathways which were previously reported to be regulated by the 
NDR (Nuclear DBF2-related) kinase COT1. Using co-immunprecipitation experiments (with MYC-tagged proteins), we have 
shown that SKB1 and COT1 physically interacted and the abundance of the 75 kDa MYC::COT1 isoform was increased in a 
Askb-1 background. On the basis of immunological detection, we propose the possible involvement of PRMTs in Arg- 
methylation of COT1. 
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Introduction 

Methylated derivatives of arginine were identified 45 years ago 
[1]. Since then, protein methylation has been recognized as a post 
translational modification that plays regulatory roles in signal 
transduction, nucleic transport, activation and repression of genes 
and mRNA splicing [2]. Arginine methylation is catalyzed by 
protein arginine methyltransferases (PRMTs) that since being 
identified [3,4] have been found in the genomes of many 
eukaryotic organisms [5]. The PRMTs are divided into four 
major types. Type I and II (conserved from yeast to human) are 
responsible for monomethylation. They catalyze asymmetric 
(aDMA) and symmetric (sDMA) dimethylarginine, respectively. 
Type III and IV PRMTs catalyze the monomethylation of 
arginine residues and the formation of monomethylarginine on 
guanidinium nitrogen, respectively. So far, Type III and IV 
PRMTs have only been found in higher eukaryotes [2] . In human 
cells there are nine genes considered to encode for PRMTs. 
Among them, only PRMT1 and PRMT3 (Type I), and PRMT5 
(Type II) are conserved through eukaryotic evolution, including 
members of the fungal kingdom [6] . 

In fungi, PRMTs have been shown to be involved in several 
cellular functions. For example, PRMT1 has been shown to be 
important for coping with oxidative stress. In Aspergillus nidulans, a 
concentration of 3 mM H 2 0 2 resulted in more than 50% growth 



retardation of ArmtA. At 6 mM H 2 0 2 conidiation was completely 
abolished [7] . PRMT 1 was also shown to be important for hyphal 
development in Coprinopsis cinerea [8] and for pathogenicity in 
Fusarium graminearum, where DON mycotoxin production and 
virulence on flowering wheat heads were reduced, and increased 
sensitivity to oxidative and membrane stresses was observed [9]. 
PRMT5 was found to play a role in cell morphology in 
Schizosaccharomyces pombe [10] and Saccharomyces cerevisiae [11], in 
the increased sensitivity to oxidative stress and elevated temper- 
ature in A. nidulans [7] and required for full virulence in F. 
graminearum [9] . So far, deletion strains of PRMT3 have not been 
shown to affect morphology in any of the fungi analyzed [7,9] . 

The NDR (Nuclear Dbf2-related) protein kinases are part of the 
AGC kinase family and belong to the Ser/Thr kinase subgroup 
[12]. NDR kinases are essential components of important cellular 
processes such as morphological switches, cell proliferation, 
mitotic exit and apoptosis in a variety of eukaryotic organisms, 
and are regulated by co-activators of the MOB (Mps 1 -one binder) 
family [13]. JV. crassa COT1 is the founding member of the NDR 
kinase family, and a single amino acid substation (within the 
catalytic domain) results in slow hyphal elongation accompanied 
by hyperbranching at temperatures above 32°G [14—16]. COT1 
kinase was shown to be regulated both transcriptionally [17] as 
well as post transcriptionally, where phosphorylation was shown to 
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be involved in the COT1 -dependent regulation of hyphal 
elongation, branching and of asexual development [15]. In S. 
pombe, the COT1 homologue Orb6 has been shown to physically 
associate with Skbl, the homologue of the human PRMTS [18]. 

In the present study we report on the involvement of PRMTs as 
mediators of post translational modifications that affect growth 
and development in N. crassa and focus on the possible cross-talk 
between arginine methylation and phosphorylation, specific to the 
interactions between PRMTs and the NDR kinase COT1. We 
characterized all of the three predicted classical PRMT-encoding 
genes in this fungus and examined their involvement in growth 
rate, branching, sexual reproduction and cell wall integrity. Our 
results indicate that COT1 undergoes methylation on arginine 
residues, and that there is a physical and as well as a genetic 
interaction between SKB1 and COT1. 

Materials and Methods 

Strains, media and growth conditions 

General procedures and media used in the handling of N. crassa 
have been previously described [19] or are available through the 
Fungal Genetics Stock Center, (www.fgsc.net/Neurospora/ 
NeurosporaProtocolGuide.html). N crassa strains used in this 
study are listed in Table 1 . Strains were grown in either liquid or 
solid (supplemented with 1.5% agar) Vogel's minimal medium 
with 1.5% (w/v) sucrose. When required, the medium was 
supplemented with 100 ug/ml Histidine (Sigma Aldrich), 10 ug/ 
ml hygromycin B (Calbiochem, Riverside, CA) or 10 ug/ml 
Nourseothricin (Werner BioAgents, Jena, Germany). For antifun- 
gal inhibition experiments, the medium was supplemented with 
either 6.5 mM Polyoxin D (Kaken Chemical Company, Tokyo, 
Japan), 4.75 |J,M Caspofungin (Merck Research Laboratories, 
Rahway), 12.23 mM Voriconazole (Pfizer Inc., Sandwich, En- 
gland) or NaCl (0.5-1.2 M) (Merck Research Laboratories, 
Rahway). 

Deletion of skb-1 was performed using the Neurospora genome 
project Gene Knockout Kit (obtained from the FGSC) according 
to Colot et al. [20]. Briefly, the construct for homologous gene 
replacement of NCU01613.3 was provided by the Neurospora 
knockout project. The cassette was PCR-amplified using the KO- 
F and KO-R primers and was transformed into the N. crassa Amus- 
52 strain, by electroporation. Transformants were screened for 
their ability to grow in medium containing hygromycin B, and 
crossed with a wild-type strain (FGSC#987). Progeny from these 
crosses were screened for their lack of ability to grow on Basta (to 
eliminate the Amus-52::barn background). The deletion mutation 
was verified by Southern blot and PCR analyses. 

A MYC-tagged version of SKB1 was constructed, utilizing the 
genomic skb-1 coding sequence along with upstream sequence that 
contained the promoter. This skb-1 fragment was obtained by 
PCR amplification with primers 5'-TCCAAGGAGACGTCAA- 
GACCG and 5 'GGTGTAAGTGGGC ATC GGAAATAA T. 
The amplicon was inserted into a pDrive vector (Qiagen, Hilden, 
Germany) to create pMEl. The promoter region of skb-1 was 
excised from pMEl by Hindi and Hindlll and inserted into 
pBluescript (Stratagene) to create pME5. The predicted skb-1 
ORF was amplified from pMEl using the primers 5'-AAT- 
GAGGGGTCACAGTGATAAGCTT and 5'-CGCGGATCCC- 
GACCTAG CCAGCCTTGC and inserted into pME5 using 
Hindlll and BamHl to create pME9. The incorporation of the 
6MYC tag resulted in the introduction of an Nsil restriction site at 
both ends of the product, which was ligated into a pDrive vector to 
create pME8. The 6MYC tag was inserted at the 5' region of skb-1 
in pME9, at the Nsil site, to create pMElO. The correct 



orientation and context of the insertion was verified by sequenc- 
ing. pMElO was linearized by BamHl and co-transformed by 
electroporation into N. crassa along with a 3.3-kb Xbal fragment 
from pMP6, containing a hygromycin-resistance cassette. The 
colonies were screened for their ability to grow on hygromycin and 
Southern and western-blot analyses (using anti-Myc antibodies) 
verified the presence of the MYC-tagged SKB1 with the expected 
molecular mass. To verify proper activity of the fused 
SKB1::MYC protein, pME9 and pMElO were co-transformed 
with pNV15 [21], containing a Nourseothricin resistance cassette 
for phenotypic complementation of Askb-1. 

For constructing COTl::GFP, a unique Ndel site was intro- 
duced into pCZ13 [which is the genomic Small EcoRI cot-1 
fragment from pOY18 [16] ligated into pUC118] upstream of 
the stop codon, by using the primers 5'-CAA CTT CCG ACA 
TAT GCG TTA C-3' and 5'-GCA TAT GTC GGA AGT TGT 
TGT C-3' to create pCZ13-NdeI. The pCZ13-NdeI, containing 
the cot-1 wild-type allele, was capable of complementing the cot-1 ts 
mutant. The sGFP sequence was PCR-amplified from pMF272 
[22] using the primers: 5'-GGT TAA TTC ATA TGG TGA 
GCA AG-3' and 5-CAT ATG GAA TTC TTA CTT GTA CAG 
CTC-3', while introducing an Ndel restriction site at both ends of 
the product, and ligated into a pDrive vector (QIAgen, Hilden, 
Germany) to create pCZ91. Both pCZ91 and pCZ13-NdeI were 
digested with Ndel and the gfp encoding fragment was ligated into 
pCZ13-NdeI to create pCZ92. Proper integration of the gfp tag at 
the 5' of cot-1 CDS was verified by sequencing the resulting 
plasmid. The resulting pCZ92 plasmid was linearized with Xhol 
and co-transformed into the cot-l;Amus-52;his-3 strain along with 
the Hindlll I BamHl genomic fragment that can complement his-3 
auxotrophy. Transformants were screened for their ability to grow 
in medium lacking L-His, at 34°C. The resulting COTl::GFP 
transformants were screened by PCR to verify homologous 
recombination. Transformants carrying a single GFP tagged cot- 
1 allele at cot-1 locus were further crossed with N crassa cot-1 ts, and 
screened for their ability to grow at 34°C without L-His and for 
their sensitivity to Basta (to eliminate the Amus-52::bar R back- 
ground). Western-blot analysis using rabbit polyclonalanti-GFP 
antibodies ab290 (Abeam, Cambridge, United Kingdom) verified 
the presence of GFP-tagged COT1 with the expected molecular 
mass. 

Determining growth rate, length between branches and 
conidiation rate 

For growth-rate measurements, race tubes containing Vogel's 
sucrose minimal medium were inoculated with 10 ju.1 of a conidial 
suspension (2 x 10 6 conidia/ml). The race tubes were incubated for 
several days at 34°C and growth was measured twice daily. 

To determine the length between branches, glass slides covered 
with Vogel's sucrose minimal medium were inoculated with 
conidia of the different N. crassa and incubated overnight at 34°C. 
The edges of the growing colonies were observed microscopically 
and the entire colony edge was documented in ca. 20 photographs. 
Each photograph was analyzed using ImageJ 1.37V (Rasband, 
W.S., U.S. National Institutes of Health, Bethesda, MD, http:// 
rsb.info.nih.gov/ij/, 1997-2006). The length between each two 
successive branches of leading hyphae was determined and at least 
400 measurements were performed for each colony. The entire 
data set of each strain was used to calculate the average length 
between branches. The presented data (and calculated standard 
error) is based on the averages of at least three independent 
experiments performed with each strain. 

To determine conidiation rates, strains were grown in a 100 ml 
Erlenmeyer flasks containing 20 ml solid Vogel's sucrose minimal 
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Table 1. Neurospora crassa strains used in this study. 



Strain 


Genotype 


Source 


Wild-type 


74-OR23-1 A 


FGSC#987 


Wild-type 


ORS-SL6 a 


FGSC#4200 


Amus-52;his-3 


Am us-52::barR;his-3 


FGSC#9720 


cot-1 ts 


cot-1 (C102t) A 


FGSC #4065 


cot-1 ts 


cot-1 (C102t) a 


FGSC #4066 


MYC::COT1 


myc(start)::cot-1 


[15] 


Amob-2a 


Amob-2a hygR 


FGSC #11296 


Amob-2b 


Amob-2b hygR 


FGSC #13575 


Aamt- 1 


Aamt-1 hygR 


FGSC #11848 


Aamt-3 


Aamt-3 hygR 


FGSC #11832 


Askb-1 


Askb-1 hygR 


This study 


Aamt-1; Aamt-3 


Aamt-1; Aamt-3 hygR 


This study 


Aamt-1; Askb-1 


Aamt-1; Askb-1 hygR 


This study 


Aamt-1 ;cor-;(MYC::COT1) 


Aamt-1; myc(start)::cot-1 hygR 


This study 


Aamt-3;cot- 1 (MYC::COT1 ) 


Aamt-3; myc(start)::cot-1 hygR 


This study 


AsM>-/;cof-7(MYC::COT1) 


Askb-1; myc(start)::cot-1 hygR 


This study 


Aamt- 1;Aamt-3;cot- 1 (MYC::COT1 ) 


Aamt-1 ; Aamt-3; myc(start)::cot-1 hygR 


This study 


Aamf- 1;Askb- 1;cot- 1 fMYC::COT1 ) 


Aamt- 1;Askb- 1; myc(start)::cot- 1 hygR 


This study 


Aamt-1 ;cot-1 ts 


Aamt-1; cot-1(C102t) hygR 


This study 


Aamt-3;cot-1 ts 


Aamt-3; cot-1(C102t) hygR 


This study 


Askb-1 ;cot-1 ts 


Askb-1; cot-1(C102t) hygR 


This study 


Aamt-1; Aamt-3;cot-1 ts 


Aamt-1; Aamt-3; cot-1(C102t) hygR 


This study 


Aamt-1; Askb-1 ;cot-1 ts 


Aamt-1; Askb-1; cot-1(C102tj hygR 


This study 


COT1::CFP 


cot-1 ::sgfp+ 


This study 


Askb-1 ;COT1::GFP 


Askb-1; cot-1 ::sgfp+ hygR 


This study 


SKB1MYC 


myc.-.skb- 1 


This study 


Askb-1; SKBIMYC 


Askb-1; myc::skb-1 NatR 


This study 



doi:1 0.1 371 /journal.pone.0080756.t001 



medium. The flasks were incubated for several days at 34°C until 
conidia were produced and matured. Conidia were collected by 
adding 50 ml of distilled water and vigorously shaking the flasks. 
The conidial suspensions were filtered through cheesecloth and 
subsequently centrifuged for 5 min at 3,000 g and the supernatant 
completely removed. The conidia were then re-suspended in 
50 ml of distilled water and the number of total CFUs was 
determined by use of a hemocytometer. 

All the results presented are the average of at least three 
independent experiments, in which minimum of triplicate counts 
were performed, with 3-5 different siblings of each strain. 

Protein extraction and immunoblotting 

Mycelial samples were frozen in liquid nitrogen, pulverized, and 
suspended in lysis buffer (1 M sorbitol, 10 mM HEPES [pH 7.5], 
5 mM EDTA, 5 mM EGTA, 5 mM NaF, 0.1 M KC1, 0.2% 
Triton X-100, complete protease inhibitor mixture [Roche 
Applied Science, Mannheim, Germany). The samples were 
homogenized by 10 strokes of pestle A in a Dounce homogenizer. 
The homogenates were centrifuged for 40 min at 10,000 g, and 
the supernatants were recovered and stored at -70°C until 
analyzed. Proteins were separated by NuPAGE® Bis-Tris 4—12% 



SDS-PAGE (Invitrogen, Carlsbad, Ca) or by SDS-PAGE 7.5% 
gels. Western blotting was performed by standard procedures [23]. 
Antibodies used throughout this study included oc-cMYC (Santa 
Cruz, CA), SYM24 and SYM10 (EMD MiUipore Corporation, 
Billerica, MA, USA), f}-Tubulin (Abeam, Cambridge, United 
Kingdom), a -COT1 [14] and a-Phosph-Ser CSK the S. cereviciae 
homolog of COT1 [24]. The secondary antibodies used were 
either goat anti-rabbit IgG-HRP (Santa Cruz, CA) for SYM24, 
SYM24, fi-Tubulin and a-Phosph-Ser CSK or goat anti-mouse 
IgG-HRP (Santa Cruz, CA) for (X-cMYC. 

Immunoprecipitation of the tagged proteins by either oc-cMYC 
or a-COTl antibodies was performed according to Maerz el al. 
[25]. For protein detection in polyacrylamide gels the Pierce Silver 
Stain Kit (Rockford, IL) was used, according to the manufacturer's 
protocol. 

Microscopy 

Light and fluorescent microscopy was performed with a Zeiss 
Axioscope microscope equipped with a Nikon DXM1200F digital 
camera. 
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Results 

Identifying PRMT-encoding genes and constructing 
deletion mutants 

The genome of JV. crassa contains three PRMT genes, which 
have typical arginine methyltransferase domains. NCU07459 
(designated amt-1), is homologous to the human PRMT1, 
NCU01669 (amt-3) is homologous to PRMT3, and NCU01613 
(skb-1) is homologous to PRMT5 [and the S. cerevisiae, as well as S. 
pombe Shkl -binding protein (Skbl)]. While knockout strains of amt- 
1 and amt-3 (the only two predicted type I arginine methyltrans- 
ferase genes in JV. crassa), were available, we produced an skb-1 (a 
predicted type II arginine methyltransferase) mutant as a part of 
this study. To do so, we used a gene replacement cassette (obtained 
from the FGSC), and transformed the JV crassa Amus-52 strain. 
Hygromycin resistant transformants were isolated and disruption 
of the skb-1 gene was confirmed by PCR and Southern blot 
analysis. All the PRMT-encoding deletion strains produced in this 
study were viable, although radial growth was slower than the wild 
type (see below).The amt-l;skb-l double mutant, also exhibited 
abnormal hyphal morphology (Fig SI). 

In order to determine whether inactivation of the PRMT- 
encoding genes conferred changes in arginine methylation, we 
performed western blot analyses, with specific anti-methyl arginine 
antibodies, on protein extracts obtained from the mutants (Fig. 1). 
Antibodies used were either SYM24, that identify asymmetric 
dimethylarginine, or SYM10 that identify symmetric dimethylar- 
ginine. Use of these different antibodies allowed us to distinguish 
between the outcomes of inactivation of the different PRMT- 
encoding genes. Differences in the asymmetric dimethylarginine 
protein profile were observed in Aamt-1 and Aamt-3 (the predicted 
type I PRMTs), while no such changes were observed in the Askb- 
1 strain (Fig. 1A). At the same time, minor differences in 
symmetric dimethylarginine were observed in Askb-1, but not in 
the other single deletion strains (Fig. IB). 

In order to further examine the involvement of the PRMTs in 
the arginine methylation pattern and to identify possible overlap- 
ping functions between the different PRMTs, we crossed deletion 
strains to produce the following double mutants: Aamt- 1 ; Aamt-3 (in 
which the two identified Type I PRMT-encoding genes were 
inactivated) and Aamt- 1 ;Askb- 1 (in which a single Type I and a 
Type II PRMT-encoding gene were inactivated). The methylation 
profile of the Aamt-1; Aamt-3 mutant detected by Western blot 
analysis indicated that deletion of these two putative Type I 
PRMTs affected not only the asymmetric arginine dimethylation, 
but also the symmetric arginine methylation profile. Furthermore, 
the Aamt- 1 ;Askb- 1 asymmetric dimethylarginine profile was similar 
to that observed in the amt-1 deletion strain, and not that observed 
in Askb-1 protein extracts. Based on the methylation profiles of the 
single and double deletion mutants tested, we concluded that 
arginine methylation was affected in all mutant combinations and 
that all three genes most likely encode PRMTs. Moreover, the 
results obtained in the western blot analyses suggest the presence 
of some reciprocal interactions between different PRMTs. 

Deletion of PRMT-encoding genes affects fungal 
morphology and physiology 

In order to further understand the role of the PRMT-encoding 
genes in the development of JV crassa, we monitored hyphal growth 
and integrity as well as asexual and sexual spore formation in the 
different mutants. 

All the PRMT deletion mutants exhibited significantly reduced 
hyphal elongation rates, though to different extents (Table 2). 
Aamt-1 grew at a rate that was approximately 30% slower than the 
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Figure 1. Changes in arginine methylation profiles in PRMT- 
encoding gene deletion strains. Western-blot analysis of total cell 
extracts probed with: (A) SYM24 antibodies used to identify asymmetric 
dimethylarginine and (B) SYM10 antibodies used to identify symmetric 
dimethylarginine. The most prominent proteins in which differences in 
the Arg-methylated signals were detected are marked by arrows. 
doi:10.1371/journal.pone.0080756.g001 

wild type, while Askb-1 and Aamt-3 exhibited a milder, albeit still 
significant, reduction (~16% and 8%, respectively) in hyphal 
elongation rates. Aamt-3 was also affected in its branching 
frequency, as distance between branches increased by 15%. 
Amongst the different types of asexual conidia produced by JV. 
crassa, under standard growth conditions, macroconida are the 
most abundant (over 90% of CFUs; [26]). Deletion of skb-1 had a 
very marked effect on conidiation, almost doubling the number of 
macro-conidia produced, while Aamt-1 and Aamt-3 had no 
observable effects on this developmental process. We examined 
conidial viability and germination rates in all the PRMT mutants 
and found no significant effect of the deletion on these parameters. 
Furthermore, all three mutants did not exhibit defects in sexual 
reproduction (as has been observed in other fungi; [7,9]). To 
further determine the possible effects conferred by deletion of the 
PRMT-encoding genes on the hyphal cell, we examined colony 
growth rates upon exposing the mutants to a variety of stresses 
known to affect cell integrity (Fig 2). We found Askb-1 and Askb- 
l;Aamt-l to be less sensitive to the chitin synthase inhibitor 
Polyoxin D (Fig. 2B), and no inhibition was detected in these 
mutants at the EC 2 5 wild type level tested. In contrast, no change 
in sensitivity was observed when the mutants were cultured in the 
presence of the glucan synthesis inhibitor Gaspofungin (data not 
shown), suggesting that this PRMT may be involved in regulation 
of cell wall biosynthetic machinery, specifically chitin synthases. 
The hypersensitivity of amt-1 (but not Aamt-3 or Askb-1) to 
Voriconazole (Fig. 2C) suggests ergosterol synthesis or membrane 
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Table 2. Morphological consequences of deletion in PRMT-encoding genes. 





Developmental process 


Wild type 


Aamt- 1 


\amt~3 


Askb-1 


Aamt-1; Aamt-3 


Aamt-J; Askb-1 


Hyphal growth rate (mm/h) 


5.3±0.07 A 


3.54±0.18 DE 


4.88 ±0.04 B 


4.47 ±0.1 3 C 


3.89±0.19 D 


3.26 ±0.07 E 


Length between branches (|im) 


162±5 B 


154±4 B 


183±5 A 


159±4 B 


154±3 B 


174±4 A 


Total CFU (*10 6 ) 


72±4 B 


68±9 B 


79±9 B 


134±7 A 


73±4 B 


127±8 A 



Growth, branching and conidiation of single and double deletion of PRMT-encoding genes. Growth rate (mm/h), length between branches (urn), conidiation, (total CFU 
*10"6). All experiments were conducted at 34°C, in triplicates. Different letters were used to mark statistically significant difference (P = 0.05). 
doi:1 0.1 371 /joumal.pone.0080756.t002 



maintenance may also be affected following the impairment of 
PRMTs. 

SKB1 physically interacts with COT1, and affects its 
protein abundance 

In JV, crassa, the NDR kinase-encoding gene cot-1 was previously 
described to be involved in hyphal elongation and branching, 
condiation and cell wall integrity [15,16,27], all of which seem to 
be affected, at various levels, by deletion of the PRMT-encoding 
genes studied here, suggesting a possible functional link between 
PRMTs and COT1. Furthermore, the presence of a physical 
interaction between Skblp and Orb6p (a COT1 homologue) were 
reported in S. pombe [18]. Thus, in order to determine if SKB1 and 
COT1 physically interact in JV. crassa, we generated an 
SKB1::MYC fusion construct, under the control of the native 
skb-1 promoter (pMElO) and inserted it into the wild-type and skb- 
1 strains. The SKB1::MYC construct successfully complemented 
the Askb-1 mutant and normal conidiation levels and growth rate 
were restored. A co-immuneprecipitation procedure was conduct- 
ed using both anti-Myc and anti-COTl antibodies, and the 
resolved proteins were subjected to western blot analysis using 



/ WT 


\amt-1 


lamt-3\ 


\lskb-1 


\amt-1; 


Aamt-1; J 




\amt-3 


\skb-lj 





Figure 2. Effect of Polyoxin D and Voriconazole on growth of N. 
crassa PRMT mutants. All strains were cultured on sorbose (2%)- 
amended (which restricts colony growth) medium at 34 r C. (A) control 
(B) with PolyoxinD (6.5 mM) (C) Voriconazole (12.23 mM). 
doi:1 0.1 371 /journal.pone.0080756.g002 



anti-Myc antibodies (Fig 3, left panel). The results obtained 
confirm the existence of a physical interaction, in vivo, between 
SKB1 and COT1. The specificity of this interaction is substan- 
tiated by the significantly lower abundance of SKB1 in a cot-1 ts 
background (Fig. 3, right panel), in which COT1 is less stable [27]. 
To further explore whether this physical interaction may affect 
COT1 intercellular localization and protein levels, we crossed the 
skb-1 deletion mutant with MYC::COTl, and COT::GFP strains. 
No changes were observed in COTl::GFP localization in the 
Askb-1 background, and COT1 was clearly observed at the hyphal 
apical region (Fig. S2), as has been previously described [28]. 

However, the deletion of skb-1 did affect the abundance of 
COT1 protein levels. Recendy, 3 isoforms of COT 1 (tagged with a 
6xMYC tag) have been identified and characterized, referring to 
protein weights of 70, 75 and 80 kDa [29]. The 70 and 80 kDa 
proteins correspond to the two previously described COT1 
isoforms [14]. In the Askb-1 background, the abundance of the 
mid-size isoform (75 kDa), was significantly elevated (Fig. 4A). We 
explored this further, and looked at COT1 protein levels during 
different developmental stages, and found that the absence of skb-1 
results in altered COT1 protein levels as early as during conidial 
germination (Fig. 4B). 

COT1 has been shown to interact with co-activators from the 
MOB family, MOB2A and MOB2 [25]. In order to determine 
whether this interaction is dependent on the physical interaction 
between COT1 and SKB1, we preformed co-immunopercipita- 
tion experiments with proteins in a Askb-1 background, and found 
that the interaction between COT1 and MOB2A or MOB2B (Fig. 
5) was not altered. These results may indicate that SKB1 is not 
required for COT-MOB2A/B binding, yet may be a part of the 
COT1-MOB2 complex. 

COT1 is methylated on Arginine residues 

Based on the detection of a physical interaction between COT1 
and SKB1, we examined whether COT1 undergoes arginine 
methylation. COT1 proteins were immunopercipited from 
MYC::COTl strain extracts and the 3 isoforms (80, 75 and 
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Figure 3. Co- Immunoprecipitation of MYC-SKB1 and COT1. 

Corresponding proteins were detected with anti-Myc, in wild-type and 
cof-7 ts backgrounds. The experiment was performed under restrictive 
temperature (36 C C). "|Cell extract" represents control without immu- 
noprecipitation. 

doi:10.1371/journal.pone.0080756.g003 
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Figure 4. Deletion of skb-1 affects COT1 protein level. Western- 
blot analysis of total cell extracts probed with anti-Myc antibodies (A) 
Single and double PRMT deletion strains in a MYC::COT1 background (B) 
MYC::COT1 in aAskb-1 background at two developmental stages, p 
tubulin levels are presented as controls for each developmental stage 
examined.. MYC::COT1 isoforms are marked by arrows. 
doi:1 0.1 371 /journal.pone.0080756.g004 

70 kDa) probed for arginine methylation using the antibodies 
described above. When the SYM24 antibodies (specific for 
asymmetric dimethylarginine residues) were used, all three 
isoforms were detected (Fig. 6A). The SYM10 antibodies reacted 
only with the 80 kDa isoform (Fig. 6B). These analyses demon- 
strate that all GOT1 isoforms are methylated on Arginine residues, 
though the nature of the methylation differs amongst them. We 
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Figure 6. Arginine methylation on MYC::COT1. Immunoprecipi- 
tation of MYC::COT1with anti-Myc antibodies probed with: (A) SYM24 
antibodies (detecting asymmetric dimethylarginine) and (B) SYM10 
antibodies (detecting symmetric dimethylarginine). MYC::COT1 isoforms 
are marked by arrows. 
doi:1 0.1 371/journal.pone.0080756.g006 

used the same approach (using the myc-cot-1 strain) to determine 
whether the Arg methylation state on the different COT1 isoforms 
is dependent on SKB1 or the other PRMT-encoding genes. 
Unexpectedly, no observable changes in COT1 arginine methyl- 
ation signals were detected in extracts obtained from single 
PRMT-encoding gene mutants (taking into consideration that an 
increase in COT1 abundance was detected in a Askb-1 
background, as shown in Fig 4). However, in the Aamt- 1 ;Aamt-3 
and Aamt- 1; Askb-1, the SYM10 antibodies could detect both the 
80 and 70 kDa isoforms of COT1 (Fig. 6B). As COT1 was 
detected with both asymmetric and symmetric dimethylarginine 
antibodies, we hypothesize that it is likely that more than one 
PRMT is involved in the process. 

In the light of these results we explored the possibility that other 
PRMTs, besides SKB1, may be involved in COT1 regulation. It 
seems that COT1 protein abundance is elevated in the Aamt- 
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Figure 5. Deletion of skb-1 does not affect COT1 interaction with MOB2 co-activators. Immunoprecipitation of MYC::COT1 with anti-Myc 
antibodies. Proteins obtained were resolved on a gradient gel and subjected to silver staining. MOB2A and MOB2B co-precipitants are marked. 
MYC::COT1 isoforms are marked by arrows. * marks the IgG. 
doi:1 0.1 371 /journal.pone.0080756.g005 
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l;Askb-l mutant, but no apparent changes occurred in the other 
deletion strains (Fig. 4A). Furthermore, utilizing specific antibod- 
ies, we have not found significant differences in the ability of 
Ser417, which is the COT1 autophosphorylation site [15,28], to 
undergo phosphorylation (Fig. S3). In addition, no changes were 
found in the ability of COT1 to interact with the MOB2 co- 
activators (Fig. 5). Taken together, it appears that if Arg 
methylation in COT1 is dependent on PRMTs, COT1 autophos- 
phorylation and MOB2A/B protein binding are not affected by 
the lack of these enzymes. 

PRMTs genetically interact with COT1 

Since we suspect that at least two PRMT are involved in COT1 
arginine methylation, we examined if skb-l or amt-1 1 ' amt-3 can 
genetically interact with cot-1. To do so we used the cot-1 ts mutant, 
in which changes are conferred in COT1 isoform abundance 
when the strain is grown at a temperature above 32°C, conditions 
in which growth arrest and hyper branching occur [14,16]. We 
crossed the single or double PRMT-encoding gene deletion strains 
into the cot-1 ts background, and examined the morphology and 
growth rate at a semi-restrictive temperature of 30°C. The 
morphology of the all single and double PRMT-encoding gene 
deletion strains in a cot-1 ts background was different from that 
observed in cot-1 ts (Fig. 7). This was mainly characterized by a 
significant increase in growth rates of all the strains at a semi- 
restrictive temperature (in comparison to permissive conditions; 
Table 3). At permissive temperatures no differences were observed 
in relation to the parental strains in a non-cot-1 ts background (data 
not shown). 

The cot-1 mutation is environmentally suppressed at restrictive 
temperature by osmotic or oxidative stresses [27]. We examined 
the ability of the various Aamt-1 /3;cot-l and Askb-l;cot-l strains to 
undergo phenotypic suppression by NaCl. In both Askb-1 and 
Aamt-3 backgrounds, the degree of phenotypic suppression of cot-1 
was significantly increased, suggesting that both these PRMT- 
encoding genes may be involved in environmental sensing by 
COT1 (Fig. 8). 

Discussion 

Different PRMT deletion mutants exhibit variations in Arg 
methylation patterns 

Eukaryotic cells extend their polypeptide diversity beyond 
constrains of the encoded amino acid sequence and folding by 
means of post translational modifications. The understanding of 
the significance of these changes for fungal growth and develop- 
ment processes has risen in the last years [30]. One such 
modification is Arg methylation, carried out by protein arginine 
methyltransferases. This family of enzymes is conserved from yeast 
to human, and regulates stability, localization and activity of 
proteins [2] . Systematic analysis of the involvement of the PRMT 
family in fungal morphology was previously explored in A. nidulans 
[7] and F. graminearum [9] . Here, we analyzed the three PRMTs of 
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Figure 7. Colony morphology of PRMT- encoding gene 
deletion strains in a cot-1 ts background. Strains were cultured 
for 72 hours at 30 C C (a semi-restrictive temperature). 
doi:10.1371/joumal.pone.0080756.g007 
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Figure 8. Effect of inactivation of Aamt-3 or Jskb-l on growth of 
cot-1 ts. cof-7 ts, grown in the presence of NaCl, exhibits a suppressed 
colonial phenotype. The colonial phenotype is further suppressed upon 
inactivation of amt-3 or skb-1. Strains were cultured for 72 hours at 34°C 
(cot-1 ts restrictive temperature) with 0.5, 0.75, 1 or 1.25 M NaCl. Colony 
area is presented as percent of control. * and + indicate a significant 
difference (p<0.05) of colony area in the Aamf-3;cof-7 or Askb-1;cot-1 
backgrounds, respectively. 
doi:10.1371/journal.pone.0080756.g008 



Table 3. Growth rate of PRMT-encoding genes deletion strains in a cof-7 ts background. 





Radial growth rate 


Aamt-1 


Aamt-3 


Askb-1 


Aamt-1; Aamt-3 


Aamt-1; Askb-1 


25 C (% from cof-7) 


77.6%±3.1% 


95.4%±2.5% 


86.4%±0.8% 


64.8% ±0.1% 


64.2%±0.4% 


30°C (% from cof-7) 


89.4%±2.1% 


128% ±4% 


105.8%±3.3% 


88.6%±4.9% 


83.9% ±2.2% 



Growth rate is presented as % of cof-7 ts. Experiments were conducted at 30 C C or 25 'C (semi-restrictive or permissive temperatures, respectively), in triplicates. 
doi:1 0.1 371 /journal.pone.0080756.t003 
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JV. crassa, with emphasis on the possible involvement of PRMTs in 
the characteristic rapid hyphal growth of this fungus [31]. 

The three JV crassa genes that contain the arginine methyltrans- 
ferase domain were identified as: NCU07459 (amt-1), NCU01669 
(amt-3), which comprise the two predicted type I PRMTs, and 
NCU01613 (skb-1), the predicted type II PRMT. Their phyloge- 
netic associations have been previously described [9]. In addition 
to the three highly conserved PRMTs, another, protein 
(NCU08111), which may be a putative type IV PRMT, but 
doesn't contain the conserved methytransferase domain, is present 
in the JV". crassa genome. In F. graminearum inactivation of the 
homologue of this gene, which is structurally significantly different 
than the cross-kingdom classes of PRMTs [5] , had no observable 
phenotypic consequences [9]. The enzymatic characteristics of this 
enzyme and its involvement in cellular functions in JV crassa and 
additional fungi have yet to be determined. 

The amenability of JV. crassa to genetic manipulations allowed us 
to construct double deletion mutants of genes from the PRMT 
gene family: Aamt-1; Aamt- 3, that is similar in nature (a double 
Type I PRMT mutant) to a mutant previously examined in F. 
graminearum [9] and Aamt- 1 ;Askb- 1 (combining a Type I and Type 
II PRMT inactivation) that to the best of our knowledge was 
produced and examined for the first time during this study. 

We used different anti methyl Arg antibodies to explore the 
profile of proteins that are methylated on arginine residues in 
PRMT deletion mutants. Our results verified the involvement of 
the predicted genes in this process. Based on these analyses we 
concluded that even though both AMT1 and AMT3 perform 
aDMA, the significant differences in the aDMA profiles observed 
between the two mutants lacking these enzymes suggest that the 
their function involves different substrates in JV. crassa. We have 
also demonstrated that, as expected, SKB1 preforms sDMA. 
Examining the Arg methylation protein profiles of Aamt- 1 ;Aamt-3 
revealed, to our surprise, a significant reduction not only in 
aDMA, but also in sDMA. Since the Arg methylation profile of 
deletion mutant in two type I PRMTs has never been examined in 
other organisms, it has yet to be determined if the results obtained 
are unique to JV. crassa. One possible explanation for the observed 
changes is the presence of reciprocal regulation of the different 
PRMTs. Specifically, that amt-1 and amt-3 are involved in 
regulation of skb-1, and vice versa. Such a form of regulation 
could function at the protein level, in a manner resembling the 
way human PRMT 1,4,6,8 can undergo self Arg-methylation [32], 
or at the transcriptional level, as seen in F. graminearum where a 
44% reduction in AMT4 (the skb-1 homologue) transcript 
abundance was detected in a Aamtl strain [9]. 

Another interesting finding was in the Arg methylation profiles 
of the Aamt-l;Askb-l strain. While in the Askb-1 mutant differences 
in sDMA were observed, in Aamt- 1; Askb-1 no such differences 
were detected. This further supports our suggestion concerning 
cooperative/reciprocal regulation of enzymes of the PRMT 
family. 

Different PRMT deletion mutants affect fungal 
morphology 

In order to improve our understanding concerning the 
involvement of PRMTs in JV crassa morphology and physiology, 
we analyzed different parameters of fungal growth and develop- 
ment. All strains with single or double deletion of genes encoding 
for PRMTs exhibited reduced growth rate. This was most 
pronounced in the Aamt-1 strain 30% ~) reduction), and resembles 
that observed in F. graminearum [9], but not in A. nidulans [7]. A 
lesser, yet significant, reduction in growth rate was also observed in 
Aamt-3 and Askb-1, that so far has not been reported in 



filamentous fungi. The double deletion mutants exhibited growth 
rates that were similar to Aamt-1. This is in agreement with that 
observed in the double deletion strains of F. graminearum [9] . These 
results raise the possibility that all the PRMTs regulate growth 
through similar cellular pathways, and Aamt-1 serves as a rate 
limiting factor in this regulatory process. 

The growth and branching of JV crassa hyphae are important for 
utilization of organic substrates [33]. We identified a new role of 
Aamt-3, showing it is involved in regulation of hyphal branching, 
as inactivation of this gene significantly increased (~15%) the 
length between branches, reducing the number of branches in the 
hyphal colony. This, along with the fact that information on 
function of PRMTs in fungi, as well as in other organisms is 
lacking [5] (, warrants further investigation on the roles AMT3 
plays. 

The skb-1 gene was the only PRMT-encoding gene that was 
clearly involved in asexual sporulation, as its deletion conferred a 
doubling in the number of macro-conidia, that represent 90% of 
all a-sexual conidia in JV. crassa [26]. However, no such 
involvement has been reported in other conidia-producing fungi 
[7,9]. Moreover, observations described here suggest that skb-1 
may be involved in chitin synthesis, an important component of 
the fungal cell wall. The hypersensitivity of Aamt-1 to Voricona- 
zole suggests it may be involved in ergosterol synthesis or other 
aspects of membrane integrity. This is in line with the growth 
reduction in the presence of SDS which was previously described 
in F. graminearum, where AMT1 was suggested to be involved in 
responses to membrane stresses [9]. Interestingly, although 
deletion mutants of PRMT1 in other fungi exhibited sensitivity 
to oxidative stress, we have not observed a similar phenotype in JV. 
crassa (data not shown). 

COT1 is methylated on Arg residues 

The PRMTs in JV crassa affect morphological aspects of fungal 
growth and development (e.g., hyphal elongation, branching, 
condiation, and chitin synthesis), which were previously reported 
to be regulated by COT1 [15,16,27]. COT1 is an NDR kinase, 
which is a subclass of the AGC class of Ser/Thr protein kinases 
[13]. Members of the NDR family, as are those of the PRMT 
family, are highly conserved, from fungi to mammals. Using 
antibodies that identify methylation on Arg residues, we found the 
presence of these modifications on COT1. While all three COT1 
isoforms can undergo aDMA, only the 80 kDa isoform exhibited 
sDMA. The most significant change in arginine methylation was 
detected in the double PRMTs deletion strains, so it seems that 
COT1 is modified by more than one enzyme from this family. 
This does not rule out minor changes in methylation which may 
occur in the single mutants but are difficult to quantify on the basis 
of the analysis performed here. 

SKB1 physically interacts with COT1, and affects its 
protein abundance 

Previous work in S. pombe identified a physical interaction 
between SKB1 and Orb6p (the COT1 homologue) [18]. We 
established that a similar interaction occurs in JV. crassa, verifying 
its existence by co-immunoprecipitation. We also determined that 
deletion of skb-1 did not affect the interaction between COT1 and 
its co-activators (MOB2A/B), suggesting that the interaction 
between these activators and COT1 is not dependent on the 
presence of SKB 1 . In order to better understand the nature of 
PRMT involvement in COT1 function, we examined the 
influence of PRMT-encoding gene inactivation on COT1 protein 
levels and found an increase in the abundance of the COT1 
75 kDa isoform in Askb-1. Interestingly, this was less pronounced 
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in the Aamt-l;Askb-l strain, though a clear signal was still 
observed. As some of the features of this double mutant were 
more similar to the Aamt-1 strain, it is possible that amt-1 has a 
more marked effect on the abundance of the 75 kDa isoform, even 
in a skb-1 background. The 75 kDa COT1 isoform is suggested to 
be the result of a yet to be identified post-translational 
modification of COT1 and its abundance is reduced following 
changes in the N terminal region (NTR) of COT1 [29]. It is 
tempting to speculate that changes in the NTR (e.g., extent of Arg 
methylation or via interaction with other proteins), as conferred by 
SKB 1 , may be involved in the changes observed in abundance of 
the 75 kDa isoform. 

PRMT-encoding genes genetically interact with cot-1 

We examined PRMTs ability to genetically interact with 
COT1, using a cot-1 temperature sensitive mutant, which exhibits 
restricted, colonial growth under restrictive temperatures. The 
results obtained suggest that all deletions in PRMT encoding genes 
affect the morphology and growth rate of a cot-1 ts mutant. 
Examination of the ability of the deletion strains to undergo 
environmental suppression at restrictive temperature, which is 
typical to the cot-1 ts mutant [14], showed that growth rate of cot-1 
was increased in both Aamt-3 and Askb-1 backgrounds, when 
grown in the presence of high salt concentrations. Thus, these 
PRMTs may be linked with environmental sensing functions 
associated with COT1. 

Conclusions 

In this study we have progressed our understanding of the 
PRMT family in regulation of morphology and physiology in JV. 
crassa, and have identified previously unknown roles for these 
genes/enzymes. Construction of double deletion PRMT -encod- 
ing gene mutants suggest there is reciprocal, direct or indirect, 
regulation of enzymes in this family. From the double mutants 
analysis it seems that there is an epistatic relationship in different 
cellular processes; the deletion of amt-1 was epistatic to amt-3 in all 
morphological characteristics. We have shown, for the first time, 
that the NDR kinase COT 1 undergoes Arg methylation and 
genetically interacts with PRMTs. These data demonstrate the 
importance of methylation as a modification involved in regulation 
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Supporting Information 

Figure SI Fungal morphology of the PRMT-encoding 
gene deletion strains. Hyphal morphology at the edge of the 
colony of wild type and PRMT-encoding gene deletion strains 
cultured overnight at 34°C. Arrow marks abnormal hyphal 
growth. Bar = 100 um. 
(TIF) 

Figure S2 Deletion of skb-1 does not affect COTl::GFP 
hyphal tip localization. Localization of COTl::GFP in wild 
type (A) and in a Askb-1 (B) backgrounds. Strains were cultured 
overnight at 34°C. 
(TIF) 

Figure S3 Deletion of PRMT-encoding genes does not 
impair Ser417 phosphorylationin COT1. Anti-Ph-Ser- 

CBK1 antibodies were used to detect the phosphorylated COT1 
isoforms, following immunoprecipitation, with anti-Myc antibod- 
ies in various PRMT-encoding gene mutant backgrounds. 
MYC::COTl isoforms are marked by arrows. 
(TIF) 
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